EXPONENTIAL CONVERGENCE TO EQUILIBRIUM FOR SUBCRITICAL SOLUTIONS OF THE 

BECKER-DORING EQUATIONS 



JOSE A. CANIZO AND BERTRAND LODS 



Abstract. We prove that any subcritical solution to the Becker-Doring equations converges exponen- 
tially fast to the unique steady state with same mass. Our convergence result is quantitative and we 
show that the rate of exponential decay is governed by the spectral gap for the linearized equation, 
for which several bounds are provided. This improves the known convergence result by Jabin & Ni- 
ethammer [16]. Our approach is based on a careful spectral analysis of the linearized Becker-Doring 
equation (which is new to our knowledge) in both a Hilbert setting and in certain weighted £^ spaces. 
This spectral analysis is then combined with uniform exponential moment bounds of solutions in order 
to obtain a convergence result for the nonlinear equation. 
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1. Introduction 

1.1. The Becker-Doring equations. The Becker-Doring equations are a model for the kinetics of 
first-order phase transitions, applicable to a wide variety of phenomena such as crystallization, va- 
por condensation, aggregation of lipids or phase separation in alloys. They give the time evolution 
of the size distribution of clusters of a certain substance through the follow^ing infinite system of 
ordinary differential equations: 

W,{t), i ^ 2, (1.1a) 

OO 

J^M^feW, (l-lb) 
where 

W,{t)■=a^cx{t)c,(t)-h^xc^+x(t) I ^ I. (1.2) 

Here the unknowns are the real functions a = Ci (t) for i ^ 1 an integer, and represent the density 
of clusters of size i at time t^O (this is, clusters composed of i individual particles) . They give the 
size distribution of clusters of the phase which is assumed to have a small total concentration inside 
a large ambient phase — be it clusters of crystals, lipids, or droplets of water forming in vapor The 
numbers a.i,bi (for i > 1) are the coagulation and fragmentation coejficients, respectively, and we 
always assume them to be strictly positive, i.e., 

Gi > for integer 1, h > for integer z > 2. (1.3) 

We set &i = for convenience in writing some of the equations. If we represent symbolically by 
{i} the species of clusters of size i, then (1.1) is based on the assumption that the only (relevant) 
chemical reactions that take place are 

for integer i ^ 1. These reactions all take place both ways; the rate at which {i} + {1} {i + 1} 
occurs is QiCia, while the rate at which {i + 1} — )• {i} + {1} occurs is bi+iCi+i. These rates are 
proportional to the concentrations of the species appearing in the left hand side, in agreement with 
the law of mass action from chemistry. Note that bi appears in the system (1.1) only for i ^ 2, 
since single particles cannot break any further The quantity Wi{t) defined in (1.2) thus represents 
the net rate of the reaction {i} + {1} ^ {i + 1}. The sum 

oo oo 

g-=Y^ jQ(0) = ici{t) for all t ^ 0. (1.4) 

1=1 1=1 

is usually called the density or mass of the solution, and is a conserved quantity of the evolution. 

The system (1.1) was originally proposed in [6], though with a variation: the monomer density 
ci was assumed constant and hence equation (1.1b) did not appear The present version, where 
the total density is conserved, was first discussed in [7] and [28], and is a widely used model in 
classical nucleation theory. We refer to these works, as well as the more recent reviews [29, 25] 
for a background on the physics and applications of the Becker-Doring equations. 

We define the detailed balance coejficients Qi recursively by 



— c,(t) = iy,„i(t) 



— c,{t) = -W^{t) 



Qi = 1, Qi+i = T-^Qt * ^ 1- 

0, + l 



(1.5) 
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An equilibrium of equation (1.1) is a constant-in-time solution (with finite mass.) For a given 
z ^ 0, the sequence := Qjz' for i > 1 is formally an equilibrium of (1.1), since all of the Wi 
vanish. Hov^rever, some of these sequences do not have finite mass (and hence are not equilibria.) 
The largest possible number Zg ^ (possibly +oo) for w^hich iQiZ^ < +oo for all ^ z < 
is called the critical monomer density (i.e., Zg is the radius of convergence of the power series with 
coefficients iQi). The quantity Zs is also called monomer saturation density, hence the subscript. 
The critical mass (or, again, saturation mass) is then defined by 

oo 

:=^iQ,z: e [0,+oo]. (1.6) 
1=1 

We emphasize that both gs and Zg are completely determined by the coefficients at, bi. It is clear 
from this that the sequences {Qiz^} are equilibria for z < Zg, and {Qizl} is also an equilibrium if 
Qs < +00. These are in fact the only finite-mass equilibria [3]. Since z n- J^i^i iQi^^ is continuous 
and strictly increasing for z < Zg (and up to z = Zg if Qg < +oo,) one sees that for any finite mass 
Q ^ Qs there is a unique equilibrium with mass g. We call a solution subcritical when it has mass 
g < gg, critical when its mass is exactly gg, and supercritical when it has mass g > gg (the two latter 
cases only make sense provided gg < +oo.) 

From the above discussion one sees that for a subcritical or critical solution there exists an 
equilibrium with the same mass as the solution; this is, there exists z ^ Zg (which must be unique) 
for which 

oo oo / oo \ 

^iQ,z*-^ic,(t) =^^c,(0) = ^? fort^O. (1.7) 

4=1 4=1 \ i=l I 

When talking about a subcritical or critical solution we will often denote by z the unique num- 
ber satisfying (1.7), and refer to it as the equilibrium monomer density. On the other hand, for 
supercritical solutions there is no equilibrium with the same mass as the solution. 

The critical density gs marks a difference in the behavior of solutions: above the critical density 
a phase transition phenomenon takes place, reflected in the fact that the excess density g — gg 
is concentrated in larger and larger clusters as time passes. On the other hand, below or at the 
critical density a stationary state of the same mass as the solution is eventually reached. Since we 
are concerned with the study of subcritical solutions we will always assume that 

Zs > (equivalently, gg > 0.) (1.8) 

A fundamental quantity related to (1.1) is the free energy H{c), defined (at least formally) for 
any sequence c = {ci)i^i by 

oo , ^ 

i7(c):=^c, (^log^-lj, (1.9) 

which decreases along solutions of (1.1) (see [3]): for a (strictly positive, suitably decaying for 
large i) solution c = c(t) = {ciit))i^i of (1.1) we have 

, oo / \ / 

a „, ,,,, ^ / Ci+i \ / ciCi Ci+1 
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for t ^ Q. Since mass is conserved, if one fixes z ^ Zs then the same identity is true of the functional 
Fz = Fz (c) given by 

OQ / \ OO OO OO 

F,{c) :=5Zc,; (log;^-l) + Y.z'Q^ = H(c) -logzY^ic. + Y.''''^-^ ^1-11^ 

this is, 

4i^.(c(t)) - -D{c{t)) for t ^ 0. (1.12) 
at 

It is readily checked that Fz{z^Qi) ~ 0, i.e., Fz vanishes at the equilibrium 

1.2. Typical coefficients. Remarkable model coefficients appearing in the theory of density-conserving 
phase transitions (see [26, 24]) are given by 

+ foralH^l, (1.13) 

for some < a ^ 1, Zg > 0, 9 > and < /i < 1. These coefficients may be derived from simple 
assumptions on the mechanism of the reactions taking place; we take particular values from [24] : 

a = 1/3, ^ = 2/3 (diffusion-limited kinetics in 3-D), 

a = 0, ^ = 1/2 (diffusion-limited kinetics in 2-D), 

a = 2/3, ^ = 2/3 (interface-reaction-limited kinetics in 3-D), 

a = 1/2, PL — 1/2 (interface-reaction-limited kinetics in 2-D). 

One obtains from (1.13) and (1.5) that in this case 

i / \ -1 

aia2 . ■ . tti-i 



( 



^263... 6. V J'-^ 

w^here the product is understood to be equal to 1 for i = 1 (so Qi = 1). One can deduce from this 
that the critical monomer density is indeed the quantity Zg appearing in (1.13), and that Qs is a 
finite positive quantity. 

Another kind of reasoning that leads to a similar set of coefficients is the following: while one 
may determine the coagulation coefficients through an understanding of the aggregation mecha- 
nism and estimate that 

a, = z" (1.15) 

for some < a ^ 1 as above, the equilibrium state of a system may be obtained from statistical 
mechanics considerations: it should be 

Ci = cxp(-ej) , 

where is the binding energy of a cluster of size i (the energy required to assemble a cluster of 
size i from i monomers, in units of Boltzmann's constant times the temperature.) Depending on 
the kind of aggregates considered, this binding energy can be estimated (at least for large i) as 

where (5 is the energy released when adding one particle to a cluster, and a is related to the surface 
tension of the aggregates. The values of /x and a for various situations are still those in (1.14). 
From this we may deduce that 

= cxp {(3{i - 1) - <j{i - lY) = zl-' cxp i-a{i - 1^) , (1.16) 
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where we define Zg := exp (— /3). By (1.5) and (1.15) we finally have 

ai^e, 6, ==Zs(i-l)"cxp(cri^-(T(i-l)^) foralH^l, (1.17) 

for some O^a^l, Zs>0 and < i^i < 1. In this case Zg is still consistent with our definition 
of the critical monomer density. Observe that for large i we have i'^ — {i — \Y ^ t^i^^^, so the 
fragmentation coefficients become roughly 

bi ^ Zs ai exp {afii'^'^) at (^Zs + -fz^) , (1-18) 

which is like (1.13) with q = Zgafi. 

Both example coefficients (1.17) and (1.13) are often used in the literature, so we always state 
explicitly the consequences of our results for them. 

1.3. Main results. The mathematical theory of the Becker-Doring equations (1.1) has been de- 
veloped in a number of papers since the first rigorous works on the topic [3, 2] . Existence and 
uniqueness of a solution for all times is ensured [3, Theorems 2.2 and 3.6] provided that 

oo 

a, = 0{i), bi^O{i), '^i^c,{0) < +00. (1.19) 

i=l 

As expected, the unique solution conserves mass (this is, (1.4) holds rigorously.) Weaker, or differ- 
ent, conditions may be imposed for this to hold, but (1.19) will be enough for our purposes in the 
present paper. Since we use some results from [16], we often also require the following: 

inf flfc ~.a>0, (1.20a) 
0< Zs< +00, (1.20b) 

(1.20c) 
(1.20d) 

We notice that the conditions in [16] are slightly different from the ones above. However, the 
differences are not essential and all arguments can be carried out with the ones here, which are 
more convenient and are satisfied by the example coefficients (1.13) and (1.15). 

The relationship between the long-time behavior of supercritical solutions and late-stage coars- 
ening theories is especially interesting, and has been studied in [27, 30, 9, 24]; in particular, in 
[24] it was rigorously shown that a suitable rescaling of supercritical solutions approximates a 
solution of the Lifshitz-Slyozov equations [18]. Asjmiptotic approximations of solutions have been 
developed in [23, 22, 11, 12]. 

For subcritical solutions (the regime which we study in the present paper) it was proved in 
[3, 2] that a stationary state with the same mass as the solution is approached as time passes. The 
rate of this convergence was studied in [16] and was found to be at least like exp(— Ct^/"^) for some 
constant C. Our purpose here is to improve this by showing that in fact the speed of convergence 
is exponential (this is, like exp(— A<) for some A > 0). We give upper and lower bounds of the 
constant A, and study its behavior as the density g of the solution approaches (from below) the 
critical density Qg. An important quantity which appears in our analysis is 



B :=sup > Q,z^ > —-- e (0,+oo]. (1.21) 



lim ■ 


Qi+l 


1 








lim 


0,1+1 


= 1. 
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Observe that B depends only on the coefficients (i ^ 1) and bi (i ^ 2). We are able to show that 
subcritical solutions of (1.1) converge exponentially to equilibrium and is a good measure of 
the speed of convergence v^rhenever B is finite: 

Theorem 1.1. Assume (1.19) and (1.20). Let c = (cj); be a nonzero subcritical solution to equation 
(1.1) with initial condition c(0) such that, for some i' > 0, 

oo 

M := ^ cxp {vi) c,(0) < +oo, (1.22) 
i=i 

and define the equilibrium monomer density z > by (1.7). Then there exist explicit T> <E (0, v) and 
Xi, > such that, for any rj G (0,17), there is a number C > which depends only on g, ?/, M and 
J2Zi e'^'z'Q^ such that 

OO 

^exp{r|i)\c^{t) - z^Q^l s^Ce^p{-X^t) forallt^O, (1.23) 

1=1 

Under these conditions B is finite, and if limi^+oo = +oo we may take A* = 1/B. 

For the particular case of the example coefficients (1.13) or (1.17), if a < 2(1 — fj,) we show 
that B can be bounded by 

Ci (log ^) "'^"^ ^B^C2 (log ^) 

for z < Zs, where Ci and C2 depend only on the coefficients (aji, {bi)i (and in particular are 
independent of z.) In the case a > 2(1 — fi) we have 

Ci S =^ C2 

for z ^ Zs and some (other) constants Ci, C2 that again depend only on the coefficients. We 
observe that if z = Zg then B is finite if and only if a ^ 2(1 — ^). 

1.4. Method of proof. Our proof is based on a study of the linearization of the Becker-Doring 
equations around the equilibrium, for which we show the existence of a spectral gap whose size 
is well estimated by 1/B (see Theorem 2.14 for details). This implies exponential convergence to 
equilibrium for the linearized system, which can be extended to the nonlinear equations by means 
of techniques developed in the literature on kinetic equations, and particularly on the Boltzmann 
equation [20, 15]. 

We observe that the improvement with respect to [16] comes from the use of a different method. 
The main tool in [16] is an inequality between the free energy (or entropy) H defined by (1.9) and 
its production rate D (see (1.10)) in the spirit of the ones available for the Boltzmann equation 
[31]. As pointed out in [16], an inequality like H ^ CD for some constant C > 0, which would 
directly imply an exponential convergence to equilibrium, is roughly analogous to a functional 
log-Sobolev inequality, which is known not to hold for a measure with an exponential tail. Since 
this is the case for the stationary solutions of the Becker-Doring equation, it is believed (though, to 
our knowledge, not proved) that this inequality does not hold in general for this equation; hence, 
the following weaker inequality (this is, weaker for small H) is proved in [16] : 
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implying a convergence like cxp(— Ci^^'^). This obstacle has a parallel in the Boltzmann equation, 
for which the corresponding inequality (known as Cercignani's conjecture) has been proved not 
to hold in general, and can be substituted by inequalities like ^ CD for e > (we refer to 

the recent review [31] for the history of the conjecture and a detailed bibliography). However, 
just as for the space homogeneous [20] and the full Boltzmann equation [15] this can be comple- 
mented by the study of the linearized equation in order to show full exponential convergence. By 
following a parallel reasoning for the Becker-Doring system we can upgrade the convergence rate 
to exponential. 

Hence, our analysis is built around a study of the linearized Becker-Doring equation, which 
is new to our knowledge. We prove here the existence of a positive spectral gap of the operator 
L, defined in Section 2 as a suitable linearization of Eq. (1.1) around the equilibrium (Qi)i^i — 

{z^Qiji'^i, in different spaces: 

(1) We provide first a spectral description of L in a Hilbert space setting. Namely, we shall 
investigate the spectral properties of the operator L in the weighted space H = £'^{Q). This 
analysis is carried out with two (complementary) techniques: on the one hand, under 
reasonable conditions on the coefficients, one can show that L is self-adjoint in V. and, 
resorting to a compactness argument, the existence of a non constructive spectral gap can 
be shown. On the other hand, using a discrete version of the weighted Hardy's inequality, 
the positivity of the spectral gap is completely characterized in terms of necessary and 
sufficient conditions on the coefficients. Moreover, and more importantly, quantitative 
estimates of this spectral gap are given. 

(2) Unfortunately, as it occurs classically for kinetic models, the Hilbert space setting which 
provides good estimates for the linearized equation is usually not suitable for the nonlinear 
equation. Thus, inspired by previous results on Navier-Stokes and Boltzmann equation 
[14, 20], we derive the spectral properties of the linearized operator in a larger weighted 
£^ space. We use for this an abstract result (see Theorem 3.1) allowing to enlarge the 
functional space in which the exponential decay of a semigroup holds. This follows recent 
techniques developed in [15], though we give a self-contained proof simplified in our 
setting. The application of this theoretical result requires some important technical efforts, 
see Theorem 3.5. 

It is worth pointing out that our techniques parallel the historical development of the study of 
the exponential decay of the homogeneous Boltzmann equation. We first show by non-constructive 
methods based on Weyl's theorem that the linearized operator L has a positive spectral gap. An 
exposition of similar techniques for the linearized homogeneous Boltzmann equation can be found 
in [8]. Explicit estimates for this spectral gap were given in [5], and similar techniques for the 
extension of the spectral gap were devised in [20] and developed in [15], and used to study the 
rate of convergence to equilibrium of the homogeneous Boltzmann equation. 

In Section 2 we carry out the plan in point (1) above, and in Section 3 we carry out point (2). 
The application to the nonlinear equation and the proofs of our main results are then given in 
Section 4. 

2. The linearized Becker-Doring equations 
2.1. The linearized operator. For the whole of Section 2 we assume the following: 
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Hypothesis 2.1. We take {ai)i^i and satisfying (1.3), define {Qi)i^i by (1.5), and assume 

(1.8). We also take < z ^ Zg and set 

Q, - Q,z\ i ^ 1. (2.1) 

We also assume that 

oo 

^:=^i2(l + a, + 5i)^Q. <+oo. (2.2) 

i=i 

We remark that, when z < Zg, the sum in (2.2) is indeed finite under "reasonable" conditions 
on the coefficients (such as (1.19)). Hence condition (2.2) is important mainly for the case z = Zg. 
Also, when z ~ Zg, condition (2.2) ensures that Qs < +oo. 

The choice of < z ^ Zg corresponds to the choice of a mass < g ^ given by (1.7). The 
unique equilibrium with mass g is precisely (Qi)i^i, given by (2.1). Notice that (1.5) implies that 

a-iQiQi = h+iQi+i, i^l. (2.3) 

Consider a solution (ci)i^i of (1.1) with mass g (this is, a subcritical or critical solution). In order 
to linearize equation (1.1) around the steady state (Qi)j we define the fluctuation h = by 

c, = (1 + h,) (2.4) 

where the components of h may have any sign. Let us carry out some formal computations, and 
leave the precise definition of the linearized operator for Section 2.2. Notice that, in order for (1.7) 
to be satisfied, it is necessary that 

oo 

iQiKii) = for all t ^ 0. (2.5) 
The weak form of (1.1) reads as follows: 

oo ^ oo 

^(^i-^Q = ^ {aiCi{t)ci{t) - bi+iCi+i{t)) - 4>i- (pi) (2.6) 

i—l i—l 

for any sequence Plugging into this the ansatz (2.4) yields, for any sequence {4>i)i 

oo ^ 

2—1 

= ^ (^a,Q,Qi{l + h,{t)) (1 + hi{t)) - 6,+iQ,+i(l + h,+i{t))^ (^<^,+i - 0, - <t>^ . (2.7) 

Using (2.3) one sees that, for any i,j ^ 1, 

a^Q^Ql{l + h^{t)) (1 + hi{t)) - 5,+iQ,+i(l + h,+i{t)) 

= a,Q,Qj ih,{t) + hi{t) - h^+i{t)) + a^Q^Qlh,{t) hi{t). 

This means that the fluctuation h{t) ~ (hi{t))i satisfies 

■^M{t) = L,{h{t)) + T,{h{t), h{t)) (2.8) 
at 
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where the linear operator L is given, in wreak form, by 

oo 

1=1 i>l 



(2.9) 



for any sequences h = {hi)i, {<j)i)i, where 

W/^ := a,Q,Qi{h, + hi) - bi+iQ,+ih,+i = a^QiQi{h, + hi - K+i). (2.10) 
Alternatively, we may write L in strong form as 

Li{h) = -^^l^t+Y.Wt^. U{h) = ^^{Wt_i~Wt) (*^2). (2.11) 
The bilinear operator r(/, g) is defined in weak form by 

oo ^ 

X!ri(/,5)Qi^i = i^^aiQiQi [hgi + fi9i) {<f>i+i - 4 - 0i) 

1=1 i>l 

for any sequences / ~ {fi)i, g ~ {gi)i and {4>i)i. Alternatively, the strong form of Ti{f,g) can be 
written as 

^ oo 

ri(/,5) = -aiQifigi --Y.a,Q, {hgi+figi) (2.12) 

and 



2 

i=l 



r.(/,g) = ^(^a,_iQ,_iQi(/,_igi + /i.g,_i) -a,Q,Qi(/,.9i +/i.g.)) * > 2. (2.13) 

Neglecting in (2.8) the quadratic term r(/i, h), one is faced with the linearized problem: 

^h{t) = L{h{t)), /i(0) = /i", (2.14) 
at 

which should be understood as the linear approximation of equation (1.1) close to the equilibrium 
(Si)i5si. Our purpose in the rest of this section is to study the operator L and its spectral properties, 
thus obtaining the asymptotic behavior of equation (2.14). 

2.2. Study of the linearized operator and proof of existence of a spectral gap. The first thing 
we need to do in order to define rigorously the operator L is to give its domain. We take expression 
(2.9) as a starting point, and we denote by ^oo the set of compactly supported sequences h — {hi)i 
(this is, the sequences for which there exists > such that hi ~ for all i ^ N). 

Definition 2.2 (The operator L). Assume Hypothesis 2.1. For a compactly supported sequence 
h = {hi)i e ioo we define L{h) by the expression (2.11) (or, equivalently (2.9) for h,(j) £ £oq). 

Notice that the only infinite sum in (2.11) converges by (2.2). (Actually, the slightly weaker 
condition J2i o-^Qi < +0° would be enough to make the definition, but we keep (2.2) for simplic- 
ity). 

One can give a more compact expression of L by direct inspection, using (2.3) repeatedly: 
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Lemma 2.3. Assume Hypothesis 2.1. For any compactly supported sequence h, 

oo 

Li{h) = -cjihi + ^ 6j7ij /or i ^ 1 (2.15) 

3=1 

where Oi are defined by 

oo 

CTi = 3aiQi +^aiQi, (Ti = a,Qi + b, fori ^2, (2.16) 

i=l 

6 J are defined by 

Si 6,2 = 22 6,1 =2&2Q2-a2QiQ2, (2.17) 

Q^^^,l = QUl,^^b^Q^-a^QlQ^ for i > 2, (2.18) 

2, 6,^-l = 2,-1 6-1,. = b,Q, for i > 2, (2.19) 
and 6,j = 0/or j + 1}. 

The numbers 6,i represent the nondiagonal entries of the infinite-dimensional matrix that de- 
fines L. One can see that the only nonzero entries of this matrix are in the diagonal and in the 
first line and column. In addition, the numbers 6,j have the important property that (2i6.j)i,j5!i 
is a sjTnmetric matrix, which suggests considering the inner product with weight Qi. There is 
another reason why this inner product appears naturally: since the nonlinear equation (1.1) has a 
Lyapunov functional (see (1.9)-(1.12)), one may look at the second-order approximation of this 
functional close to the equilibrium (2i)i5si, which happens to be 



oo 



2 

1=1 

It becomes then natural to study the linearized operator L in the Hilbert space H := £'^{Q) defined 
as 

n = [h = {hi),^i I \\h\\H<^] 

with inner product denoted by (•, •). As remarked above, L becomes a symmetric operator in %. 
Another crucial property of L in this Hilbert space is that it is dissipative; this is, 

{Lh,h) ^ for all h e ioo, 

as can be readily seen from (2.9). Though L is in general not continuous on Ti, it is easy to give a 
dense subspace of "H in which it is bounded: 

Lemma 2.4. Assume Hypothesis 2.1, and also that 

Q, 

> Zs asi^ +00. (2.21) 

Qi+l 

Then, there exists C > (depending only on z and the coefficients (ai)i^\, (hi)i^\) such that 

\\L{h)\\n^C\\h\\n. (2.22) 
for any compactly supported sequence h, where 

oo 

H2 = ^'(2(1 + a^)) := [lie'H \ \\h\W^ ^ Q,{1 + a1)h} < +00} (2.23) 

1=1 

and {ai)i was defined in (2.16). 
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Remark 2.5. About condition (2.21), notice that Zg is the radius of convergence of the power series 
with coefficients iQi; hence, if the Hmit of Qi/Qi+i exists, it can only be z^. We emphasize that all 
of these conditions are met by the example coefficients (1.13) and (1.17). 

Proof. Instead of (2.22) we will show, equivalently, that 

{Lh,4>) i^CWhWuMWu forall</)eH. (2.24) 

In this proof C denotes any positive quantity depending only on z and the coefficients (ai)i^i, 
{bi)i^i, possibly changing from line to line. From (2.9) we have, using Cauchy-Schwarz's inequal- 
ity, " 



{Lh,(j)) = '^a^QiQiijii + hi~ /ij+i) ((^i+i 



/ oo \ 1/2 / oo \ 1/2 

^Ql[J2a^^Q^{h^ + h,-h,+,Yj ( ^ - 0, - 0i) M . (2.25) 

The term inside the first parentheses is bounded by 



oo oo 



O 

where we used (2.3) and the fact that Qt+i/Qi is bounded uniformly in i thanks to (1.20c). A 
similar reasoning shows that the second parenthesis in (2.25) is bounded by 

oo 

^ Qi('/'i+i - 0i - 0i) ^C\\(t)\\li, 
1=1 

which shows (2.24). □ 



The operator L is defined only on compactly supported sequences h. We will now extend it to a 
larger domain, so that it becomes a closed operator, and study its spectrum. We will consider two 
ways to do this, later shown to lead to the same result: one of them is to look at the closure L of 
L, and the other one is to look at an associated quadratic form, which will lead to an extension 
denoted by C. Let us first describe the latter. 

Expression (2.9) suggests the introduction of a sjmimetric form £ mTihy setting 

^{h,g) = -{h,Lg) = ^aiQiQi(/ij+i ~hi~ hi){gi+i - gi - gi) for h,g ^ £00, 
which can be naturally extended to the domain 



= {he'H; S{h,h) < 00} 
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by setting 

<^{h, g) = a^Q.Qi {h,+i - h, ~ hi) {g,+i - g, - ,91) for h, g £ (2.26) 
Then, one has the following 

Proposition 2.6. Assume Hypothesis 2.1. Define the space Hi by 

■Hi = {hen-, \\h\\l^^ := ^(1 + a,)Q,h^ < 00 |, 

i=i ^ 

where we recall that ct,; was defined in (2.16). 

The form S with domain Si{S) defined hy (2.26) is a closed symmetric form on %. Its domain 
3i{S) contains Hi and there exists C > (depending only on the coefficients (0^)^^! and {bi)i^2) such 
that 

\<^ih,g)\^C\\h\\nM\n, for all h,g e Hi. (2.27) 

Proof The fact that S" is symmetric and bilinear is clear Moreover, one proves as in [13, Example 
1.2.4] that Si{S)) is closed thanks to Fatou's Lemma. Let us prove (2.27), which in turn shows 
that Hi C !^{S). By a usual depolarization argument it is enough to show (2.27) for g = h. We 
have 

00 00 
<^{h,h) ^Y,a,Q,Qi{h, + hi ^ h,+iY ^ 3^a,Q,Qi(/i2 + /i2 + /i2^^) 

1=1 i=l 

00 00 

= SAhj + S^aiSiQift? + 3Y,b^Q^h^^ ^ C\\hf^^, 

i=l i=2 

where A is the quantity in (2.2). This proves the desired bound. □ 



Due to the previous Proposition, according to [13, Theorem 1.3.1 & Corollary 1.3.1], there 
exists a unique non-positive definite self-adjoint operator C onH such that C ^((f) and 

<^{h,g) = -{Ch,g) for all h e ^{C), g £ 'S{S). (2.28) 

More precisely, Si{S) = —C). It is clear that the linear operator (£, !^{C)) extends the above 
linear operator L, defined on £00, and that £00 is a core for C. It is also easy to see, under the 
conditions of Lemma 2.4, that the domain of C must include the space H2 (since £ is a closed 
operator that extends L), and that the expression of C in H2 is still given by (2.11) (since each 
of the sums converges absolutely in this space, as deduced from the proof of Lemma 2.4) or 
alternatively by (2.9) (for any e ioo). It is also easy to see that is an eigenvalue of C with 
explicit eigen-space: 

Lemma 2.7. Assume Hypothesis 2. 1 and (1.20c). Then is an eigenvalue of C, with a one-dimensional 
associated eigenspace spanned by the sequence defined by hi ~ ifor i ^ 1. 

Proof Under the condition (2.2) one sees that h = {i)i^i e H2 '!= 'S>{C). It is clear then from 
(2.28) and (2.26) that {Ch, g) = for all g e so Ch = 0. On the other hand, if there is any 

h = {hi)t^i e !^{C) such that C{h) = 0, then S'{h, h) = and consequently hi+i = h^ + hi for all 
i ^ 1. This implies that hi = ihi for i ^ 1, so /i is a multiple of h. □ 
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It is well-known [17, 13] that the Dirichlet form ^((f)) (or, equivalently, the operator 
(£, &{C))), generates a strongly continuous semigroup of contractions {St, t ^ 0} inV. — i^{Q). 
In particular, for any ~ (/i? G ^ the linearized problem 

= C,ih{t)), h,{0) = hi (2.29) 

admits a unique mild solution h{t) = {hi{t))i^i given by h{t) = Sth° satisfying 

s$ for all OO. 

Moreover, StiH) C ^((?) for any i > and, if h" = € ^(>C) then, for any t ^ 0, 

supggjo t) ||/i(s)||®(^) < oo (we refer to [13, Chapter 1] for details). In fact, for /i" G !^{C) we 
have that h is in time, h{t) g ^(-C) for alH > and (2.29) is satisfied pointwise in time. Hence, 
for /i" G ^(>C) we have 

-^JHt)f - (/i^-C/i) - - - ^a.Q.Qi(/i.+i foralH ^ 0. 

The advantage of the approach involving the quadratic form is that we obtain, in a simple way, 
a self-adjoint extension of L on which we have a lot of information. However, it is not easy to give 
explicitly the domain of the operator C. We consider also a different approach based on the 
closure of L. Since L is symmetric, it is closable and we define L as its closure: 

Definition 2.8 (The operator L). Assume Hypothesis 2.1. We define the linear operator L, with 
domain ^(L), as the closure in V. of the linear operator L (with domain £oo)- 

The following proposition gives a certain structure to L that will allow us to identify its domain, 
prove it has a spectral gap (in a non-constructive way, using Weyl's theorem) and show that in fact 
L = £: 

Proposition 2.9. Assume Hypothesis 2.1 and conditions (1.20a), (1.20c) and (1.20d). Assume also 
that z < Zg. Then the linear operator L is self-adjoint with domain ^(L) = 7^2 (as defined in (2.23)) 
and, given S > 0, it can be written as 

L{h) = L^'{h) + L^{h) for all h e ^(L) (2.30) 

such that 

(1 ) L'^ is a compact operator on %. 

(2) is a self adjoint operator with domain 'H2, and with spectrum 



where 

Am o: I 1 - 2—^ I - (5, £ := Zs/z 




and 



a := inf ak, 

k>l 



which is strictly positive due to (1.2Ga). 
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Proof. Take an integer ^ 1, to be fixed later. We define, for all i ^ 1, 

oo 

^iih) ■■=^X{min{i,j}^N} ^.i.jhj for h e -H, (2.31) 

OO 

:= -cj,h, + ^X{m\n{t.j]>N} ii.]hj for h S ^2, (2.32) 

i=i 

We recall that the notation ai and was defined in (2.16)-(2.19). The notation X{--- } represents 
a function which is equal to 1 when the condition in the brackets is satisfied, otherwise. In a 
more explicit way, 

[Lfih) =£,ishi +bihi^iX{t^N+i} + aiQihi+iX{i^N} fori > 2; 

while 

(/i) = -cTj /i, + fei X{i>Ar+i} + fli Qi ^i+i X{i>iV} (* ^ !)• (2.34) 
With calculations similar to those in the proof of Lemma 2.4 one sees that the sums in (2.31)-(2.32) 
converge, that iP : U ^ U is bounded and that there exists C > such that \\V^h\\H ^ C\\h\\H:, 
for any h <E 'H2- Moreover, both and L''^ are easily seen to be symmetric operators. One sees 
that : H — > "H is a finite-rank operator, since except for its first component Lf it only depends 
on a finite number of components of h. In particular, it is a compact operator Since it is also 
symmetric, is thus self-adjoint. 

Let us investigate now the remaining part L^^. We will prove that, for iV > 1 large enough, 
is self-adjoint thanks to Kato-Rellich's theorem [17, Theorem 4.3, p. 287]. To do so, define the 
multiplication operator T with domain ^2 by 

T,{h) = -aih, for e ^2, i ^ 1. 

Being a multiplication operator, T is self-adjoint in Ti. Define the operator S, also with domain 
"^2, as the remaining part in (2.34): 

Si{h) = hi hi^iX{i>N+i} + ai Qi hi+iX{i>N} for heH2, « ^ 1, 

so that L^^ = T + S. Clearly, S is sjnnmetric. Let us now prove that S is T-bounded with relative 
bound smaller than one. To do so, we compute 

while 

\i=l / \i>N+l / \i>N J 

\i>N / \i>N+l / 

= fE/^'2^/iA + f E '^'SA") (2.35) 

\i>N ) \i>N^\ ) 
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where 

2 ^2 2j-l 2/-i2 Qi ,2 Qi+1 



3.-1 



Observe that, since <Ti = OiQi + bi = Qi(a.i + ai_i ^g.M , and caUing £ Zg/z, 



as i — > +CX3, 



due to (1.5), (1.20c) (which implies Qt/Qt+i £ as i ^ +00) and (1.20d). Similarly, 



as I 



{i+£r 

Hence, from (2.35) we see that for any e > we can find > such that 

2V£ ?^ .^V^^ , ^ 



|S(/i)|l«^(l + 6)^(^a2Q,/,2) ^(i + e)^||T(/i)!l«^:0|lT(/i)|l„. (2.36) 



Since we are assuming z < Zg, we have £ > 1 and it is possible to choose e > such that 9 = 
2(1 + e)V£/{l + £)<!. In other words, S is T-bounded with relative bound 9 strictly smaller 
than 1. According to Kato-Rellich's Theorem L'''' = S + T with domain ^(L^) = &{T) = 712 is 
self-adjoint. 

Let us now show that the spectrum of L'^ is to the left of Am- Since a_ = inf^^i ai > 0, one 
has {T{h),h) ^ for any h e ^(T). With the terminology of [17], this exactly means that 

the multiplication operator T is bounded from above with upper bound 7t = —g_. Combining 
inequality (2.36) with [17, Theorem 4.11, p. 291], L'^ is bounded from above with upper bound 
7lm such that 7lm ^ 7t — ^'IttI = —(1 — S)a, i.e. 

(L^h, h) ^ -(1 - 9)a\\hf for all h e 9{T). 

This proves in particular point (2) of the Proposition with Am = 1 — ^' > 0. Notice that we can take 
e > arbitrarily small, so Xm can be arbitrarily close to 2^/1/(1 + £). 

Now the identity L{h) = ^^{h) + L^'(/i) is obviously true for h G £oo- Then, it is clear from 
the bounds above that the closure of L has domain "^2 and is equal to L*^ + L'^^, which shows 
(2.30). Since L is a compact self-adjoint perturbation of L'^, Weyl's Theorem ensures that L is 
self-adjoint. □ 

Corollary 2.10. Assume the conditions of Proposition 2.9. There exists Ao > (defined in a non- 
constructive way) such that 

00 

(L/i,/i) -Aol|/i|p for all he'H2 such that ^ih,Qi=Q. (2.37) 

i=l 

In addition, the operators L and C are equal. Finally, L is the generator of a Co-semigroup {St)t^o 
H such that 

00 

\\Sth\\n'^e-^"'\\h\\n eU such that ^ih,Q, = Q, t > 0. 
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Proof. Since L*-^ is a compact perturbation of L'^, Weyl's Theorem ensures that the operator L is 
self-adjoint and its essential spectrum ©oss(L) coincides with that of L'^^ in particular ©css(L) ^ 
{~oo, —Am]- Consequently, the part of ©(L) contained in (—Am, +oo) is a set of isolated eigenval- 
ues of finite multiplicity. As we also know that L is non-positive, we have & (L) C {—oo, 0], having 
only a finite number of eigenvalues in (—Am, 0]. From Lemma 2.7 we know that is an eigenvalue 
of L with eigenspace spanned by Since Am > 0, L must have a strictly positive spectral gap 

Ao, which gives (2.37). 

Since L, the closure of L, is self-adjoint, the operator L is by definition essentially self-adjoint. 
As such, it only has one possible self-adjoint extension. Since C is another such extension, we have 
L = £. The final part of the result is a classical consequence of (2.37). □ 

Proposition 2.9 does not apply in the case z ~ Zg. In order to include that case we need to 
obtain more delicate estimates: 

Lemma 2.11. Assume the hypotheses of Proposition 2.9, but take z = Zg. Define then 

Ok ~ 1 = 1 for k ^ 1. 

OfcQi Uk Qk 

and assume that 

liminf 5fc Y^Ofe > 0. (2.38) 

Then the conclusions of Proposition 2.9 still hold true for L, for some Am > 0. As a consequence, the 
conclusions of Corollary 2.10 are true (i.e., L has a positive spectral gap and L = C). 

Proof The proof follows the same steps as that of Proposition 2.9 and is based upon the splitting 
of L as L = L*^ + L'^ with L*^ and L"^ defined by (2.33) and (2.34) respectively It will consist 
in proving that, under the supplementary assumption (2.38), one can choose > 1 large enough 
such that 

(L^(/i),/i) < -Am||/i||w foralUe-Hz (2.39) 

holds true for some positive Am > 0. Once we have this, the rest of the proofs of Proposition 2.9 
and Corollary 2.10 are still valid. Indeed, arguing as in the proof of Proposition 2.9, this implies 
the existence of a positive spectral gap (defined in a non-constructive way) for L thanks to Weyl's 
theorem. One computes, for some A^ > 1 to be fixed later, 

oo 

{L^{h),h) = -Y,<^^Q^h^ + b^Q^h,h,^l + Y,a^QlQ^h^h^+l 

i=l i>N+l i>N 

oo 

= -Y,<J^Q^h^ + 2Y,<^^Ql Q^h^ h,+i (2.40) 

i=l i>N 

where we used (1.5). For each i ^ 1 we take ri > 1 (also to be fixed later), and use Young's 
inequality to deduce that 

i>N i>N+l i>N+l 

i>N+l i>N+2 * 
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again by (1.5). Since 6j = (aiQi)(l + Si) and cr,; = + a^Qi = (a,;Qi)(2 + (5.^) for any i > 1, one 
gets 



i=l j>Ar+l 



i=l i>N+l ^ 

Take N large enough so that for Si > for all i > + 1. We make the choice 

r, = l + 5,/2 fori>iV + l, 

to obtain 

N+l „2 

According to (2.38), 



w+i ^ JV+l „ 



so choosing N large enough, 

w+i 

1=1 i>N+l i—1 i>N+l 

which yields (2.39) with Am = min{a:, f } > 0. □ 
Remark 2.12. For the typical coefficients given by (1.13), with z = Zs, one sees that 

Since au = with a > 0, one sees that (2.38) holds true if and only if a > 2(1— /i). In other words, 
L has a positive spectral gap whenever a ^ 2(1 — ^). We will prove this later on in a different way, 
by a constructive argument related to Hardy's inequality which gives explicit estimates of the size 
of the gap. We will also show that, for a < 2(1 — ^), C does not have a positive spectral gap (see 
Lemma 2.22). 

2.3. Explicit spectral gap estimates in £'^{Q). Let us now study conditions ensuring that £ admits 
a positive spectral gap in the space "H ~ £^{Q)- Since £ is a self-adjoint operator in this space, it 
having a spectral gap of size Ao is equivalent to the functional dissipativity inequality 

oo 

{h,Ch) si-Xo\\h\\^ for all /i e 4o such that ^ i Q,/i, = (2.41) 

this is, 

oo oo 

AoE 2'^^? E«'2.Qi(/i.+i - h, - hif (2.42) 

i=l 1=1 

for all h e £oo such that X^i^i * Qi^i = 0. Notice that it is enough to have the inequality for 
compactly supported sequences h, since ^oo is a core for C. On the other hand, if (2.42) holds for 
all h G ^00 then it must actually hold for all sequences h, with the understanding that either or 
both sides of the inequality may be infinite. 
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We already saw in Proposition 2.9 (see also Lemma 2.11) sufficient conditions ensuring the 
existence of a positive spectral gap. However, this existence result was based on a compactness 
argument and, consequently, does not provide any quantitative information about the size spectral 
gap. In this section we adopt a different viewpoint by studjdng the inequality (2.41) directly. This 
will result in useful estimates on the spectral gap, explicit in terms of the coefficients {ak)k^i and 
(&fe)fc>i- 

Definition 2.13. Assume Hj^othesis 2.1. We call Aq the size of the spectral gap in H of the 
operator C; equivalently, Ao is the largest nonnegative constant such that (2.41) (or equivalently, 
(2.42)) holds. Notice that Ao = whenever C has no spectral gap. 

Our main estimate of Aq is the following: 

Theorem 2.14 (Bound of the spectral gap for the linearized Becker-Doring operator). Assume 
Hypothesis 2.1. Recall the definition of B given in (1.21); 

B = 

It holds that 

^ (2.43) 
which should be understood as saying that B = +oo if Aq = 0. If we additionally assume that 

oo / oo \ 

"ft^^E;^ E.C, <+oo (2.44, 
1=1 \J=i+l / 

and call 

oo 

Mai^^i'Q. (2.45) 

1=1 

(which is finite due to (2.2) j, then we also have an upper bound o/Aq." 

which should be understood as saying that B = +oo if and only if \q = 0. 

Using this it is easy to give quite general conditions on the coefficients Oj, bi such that the 
spectral gap Ao is strictly positive whenever z < Zg. In particular, the following lemma applies 
to the "tj^ical" coefficients (1.13) and (1.17). Notice that the following result has already been 
obtained in Prop. 2.9 by a completely different argument; however, the following Corollary is 
constructive, reljdng on the above Theorem: 

Corollary 2.15. Assume Hypothesis 2.1 and conditions (1.20a), (1.20c) and (1.20d). Then for any 
< z < Zg, the spectral gap Ao of CinH is strictly positive. 

Proof Due to (2.43) it is enough to show that B < +oo. Let us call Qk := z'^Qk as usual, and 

oo ^1 

TOfe := ^ Qj, ?^fe := ^ — ^ for all fc > 1. 

j=k+i j=i 



We will show that 
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^,'=-+-(^-1)- Q, (2.46) 
n.'^^^r^-lY'^, (2.47) 



which then implies that 



mfcTT-fc = —^{ukakQk) ^ -^{nkakQk) < C for all fc ^ 1 
flfc y.k £ y.k 

for some C > 0, due to the lower bound on (1.20a) and the limits (2.46) and (2.47). This implies 
that B ~ sup^,^j(TOfe"fe) < C, so we just have to prove (2.46) and (2.47). 
We prove them by using the Stolz-Cesaro theorem. Notice that 

mfe+i - TOfc _ Qk+i _ f Qk A ^ 



Qk+i - Qk Qk+i - Qk \ Qk+i J V 2 

due to (1.20c). Since is strictly decreasing and tends to as fc — > +oo, this implies (2.46). 
In a similar way, due to (1.20c) and (1.20d), 



nk+i - Uk (afe+iQfe+i) ^ _ L_ ak+iQk+i \ ^ ^ (i_ ^ ^ 



(afc+iQfe+i) ^ - {akQk) ^ {ak+iQk+i) ^ - {akQk) ^ V «fe2 
which, since Uk is strictly increasing and unbounded, implies (2.47). □ 

Let us begin the proof of Theorem 2.14. We will prove it through an inequality which is stronger 
than (2.41), and equivalent in some cases, as we will show immediately afterwards: 

Lemma 2.16. Assume Hypothesis 2.1. Define Ai as the largest nonnegative constant such that 

oo oo 

Ai^Q.(/i, -i/ii)' ^^a»Q«Qi(/i»+i -/i, -/ii)' for all h = {h,),^i. (2.48) 

i=l i=l 

Then 

— ^M^—, (2.49) 
which again should be understood as saying that Ai = if and only if B ~ +oo. 
Proof We call 

Ml 2i+ii J^i ^ a,Qi fori^l. 

Hardy's inequality from Appendix A says that B is finite if and only if the following inequality holds 
for some Ai > 0: 

/ \ 2 

OO / * \ ^ 

1=1 \]=1 J t=l 

In addition, if any of these conditions hold, we have (2.49). Hence it is enough to show that (2.48) 
is equivalent to (2.50). 

For one implication, assume that (2.50) holds. If we take any sequence h = {hi)i^i and call 

fi ■= hi+i - /ij - hi for i ^ 1 
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then (2.50) is just (2.48): notice that the term on the right hand side of (2.50) is 

oo oo 

= II 2^ {h,+i -h,- hif (2.51) 

i=l 1=1 

and the term to the left of (2.50) is 

OO / ^ \ ^ ( ^ 

i=\ \j=\ j i=\ \j=\ 

OO 

= \xY.^i^ - (i + 1) /ii)' = Ai ^ (/i, - i hif , (2.52) 

i=l i=l 

since the first term in this sum is equal to 0. 

For the reverse implication, assume (2.48). Then for any sequence / = {fi)i^i, define recur- 
sively 

hi = 0, hi := fi + hi^i + hi for i^2. 
Since (2.51) and (2.52) still hold, (2.48) imphes (2.50), and we have finished the proof. □ 

Lemma 2.17. Assume Hypothesis 2.1. There is the following relationship between inequality (2.48) 
and the spectral gap inequality (2.41): 

where for the right inequality we need to assume additionally that both is finite. (We recall that 
A/2 and M3 were defined in (2.44)-(2.45);. 

Proof Denote the right-hand side of (2.48) by D{h) for convenience (it is equal to S'Qi, h) = 
— {Ch,h) for any h G Si{C), but we rename it to include any sequence /i). Let us show the first 
inequality. For any sequence h = {hi)i^i e ^00 with J2i iQi^i = we have 



i—l i—1 i—1 i—1 i—1 

since the mixed term obtained when expanding the square vanishes due to the orthogonaHty con- 
dition on h. We have, using (2.48), 

oo 

MMn ^ Ai ^ {h, -ihif ^ D{h). 

This shows the first inequality. 

For the second inequality we assume that Ao > (since otherwise it is a trivial statement) . Take 
any sequence h = G ^oo- Calling 



oo oo 

Y.Q^{h,-^hlf =J2Q^{h^-^h)^ ^ ^2.54) 



4=1 ^ i=l 

we have 
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since the mixed term obtained from the square is zero due to the definition of h. If we call 
gi := hi — i h we see that J2iLi iQi9i — and we may use (2.41) to obtain 



(2.55) 



This deals with the first term in (2.54). For the second term, use that 



M2 \h- hi\ 



^iQiihi - ihi) 



i=l 



^iQ,;^(/lj + l - hj - hi] 

i=l j=l 



^^'^iQi\hj+i - hj - hi\ ^^\hj+i - h-i - hi\ ^ iQ^ 

i=l j=l j = l i=i+l 

00 I 00 



1/2 



< I ^ aiQi{hj+i ~ hj - hi) 



A^/^Dih)'/^ 



to obtain 



y^Q^{^h^^hlf = {h- hifAh ^ TT^ih) 



(2.56) 



Using (2.55) and (2.56) in (2.54) we finally obtain 



Y.Q^{h^-^hlf^(^ + ^)D{h) 



for every h £ £00 • It is easy to see that this implies the inequality for all sequences h, hence proving 
the second inequality in (2.53). □ 

2.4. Estimates of the spectral gap size near the critical density. In many cases of interest one 
can estimate the quantity B as z approaches Zg. This implies an estimate on the size of the spectral 
gap for the linearized Becker-Doring equations, as stated in Theorem 2.14. We remark that our 
main estimate for B, Lemma 2.22 below, applies to the coefficients (1.13) as well as (1.17). 
Let us set the notation. We call 

5fc := -log Qfc - /clog Zs forfc^l, w := log — 

z 

so that 

Qk Qkz'' = cxp {-gk - kw) for k'^l. 

Lemma 2.18. Assume that there exists fco > such that: 

(1) The sequence {gk)k is increasing for k ^ ko. 

(2) The sequence {gk/k)k is decreasing for k ^ ko. 
Assume also that: 

(3) There are constants C > and < /i < 1 such that 

gk+i-gk-Ck^^'- (2.57) 

(4) There are constants C,a > such that 

ttk ^ Ck"^ ask^ +00. (2.58) 
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Then for some numbers Ci < C2 depending only on the coefficients {ak)k and {bk)k, we have that for 
all Zs/2 < z < Zs, 



; T ^ > Q i < ; T for fc ^ 1. (2.59) 



and 

C\ ^ 1 C2 

k^OjwTk^ ^ ^ ajQj ^ fc"Qfc (w + fcA*-!) 

Remark 2.19. We give the above estimates only for z > Zs/2 for simplicity, and since we are 
interested in the behavior for z close to Zg. They can easily be derived for z > S, for any 6 > 
(with constants Ci, C2 depending on 6), but are very poor estimates for small z. 

Remark 2.20. Both example coefficients (1.13) and (1.17) satisfy all hypotheses of Lemma 2.18. 
Notice that for (1.17) the sequence gk is 

5A: = - log Zs + ^(fc — 1)^ for fc ^ 1 

while for (1.13) we have 

5fe = a log fc - log Zg + Y^ log n + j for fc ^ 1. 

j=2 \ J / 

From this it is easy to check that all the requirements of Lemma 2.18 are met. 

Remark 2.21. Notice that, for any sequence {gk)k satisfying the assumptions of Lemma 2.18, the 
following hold: 

(i) There is a constant C > such that 

gk/k < Cfc^^i for fc ^ 1. (2.61) 

This is a simple consequence of (2.57). In particular, gk/k tends to as fc — > +00. 

(ii) There exists C > such that 

gk ^ {C/^i)k^' as fc ^ +00. (2.62) 

This is easily deduced from (2.57) by using the Stolz-Cesaro Theorem. In particular, gk 
diverges to +00. 

(iii) Points (1), (2) together with (2.61) and (2.62), show that there exists fci > such that, 
for all fc ^ fci, 

gk > gj for all j fc, (2.63) 

^ ^ for all j 5; fc. (2.64) 

fc j 

All this properties will be used repeatedly in the proof of Lemma 2.18. 

Proof of Lemma 2.18. We will use in this proof a simple bound for the following geometric sum: 

for any ^ > 

exp(-fcO ^ ^ exp(-(fc - m fo^^^o^ (2 653 

Also, we will denote by C any constant that only depends on the coefficients (ak)k and {bk)k, 
and which may change from one line to the next. We will show (2.59) and (2.60) for fc ^ fci. 



EXPONENTIAL CONVERGENCE TO EQUILIBMUM 23 

Note that this is enough, since for k bounded there are always constants Ci, C2 which satisfy 
(2.59)-(2.60). 

Step 1: Proof of the upper hound in (2.59). We use two different estimates to prove this upper 
bound. The first one is 

EQj = y^cxp(-.gj - jw) sC cxp(-gA;) V'cxp(-A:w) ^ Qk C— , (2.66) 

j = k j=k j=k 

using that gk is increasing, and also Eq. (2.65). The second one is 

00 00 00 

^Qj^Y. ^M-9i - jw) ^ eM-kw) exp(-3j) CQfefc^"'^, (2.67) 

j=k j=k j=k 

due to Lemma B. 3 (see also the Remark B. 1 immediately afterwards). We deduce from (2.66) 
and (2.67) that X^jlfc Qj bounded by the harmonic mean of both right hand sides, this is, 

VQ, < (2.68) 

j=k 

Step 2: Proof of the lower bound in (2.59). 
Using that {gk/k)k is decreasing we have 



j=k j=k 



(2.69) 



w + gk/k w + kt^ 1 

j=k 



Here we have also used (2.61). 



Step 3: Proof of the upper hound in (2.60). 

We use two different bounds, in a similar way as in Step 1. The first one is 

Y QJ^^-J^ =Y^^P^9j +3w)j'°' exp(.gfc)^exp(ju;)j-" 

where we have used (2.63). Now, since exp{jw) ^ exp{{y + l)w) and ^ j~" for any y e 
(j - l,j), one has 

^cxp(ju;)j^" s; ^ / cxp((y + l)w)y-"dy = / cxp{{y + l)w)y~°' dy. 

Using then Lemma B. 2 in Appendix B (with fi — 1), we get J2j=i exp(jw)j~" ^ C exp{kw)k~°' 
from which we deduce that 
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The second bound is 

k k 

^exp(5, ^ exp(fcu;)^exp(.g,).r" s$ CQ-'k'-''-^, 

where we used Lemma B. 3 again. Taking the harmonic mean of these two bounds gives the 
desired one as before. 

Step 4: Proof of the lower bound in (2.60). 
Similarly to (2.69) we have 



using again (2.61) and (2.64). □ 

As a consequence of the previous lemma we have an estimate of the size of the quantity B in 
(1.21), useful mainly because it gives its blow-up rate as z approaches Zg, making more precise the 
above Lemma 2.11: 

Lemma 2.22. Assume the conditions in Lemma 2.18. If a > 2(1 — ^) there are numbers Ci, C2 > 
such that 

-2+-r§T7 



B ^ C2 ( log — 

z 



for z < Zg. In the case a ^ 2(1 — /i) the quantity B is hounded uniformly up to z = Zg: there exist 
numbers Ci , C2 > such that 

Ci B < C2 

for z ^ Zg. 

Proof Using the bounds in Lemma 2.18 and the fact that Qk+i/Qk is bounded from below and 
above uniformly with respect to w, we have that, for some numbers Ci, C2, 

Ci sup ^ ^ B C2 sup ^ . (2.70) 

In the case a < 2(1 — /i), the minimum of the function a; x°'{w + x''^^)^ occurs at cc = 
Cq-.^w^^/*^^^^, for Cq_p some quantity depending only on a and ^. Hence we obtain 

Ci«;"2+T^ <; B C2w~^+T^, 

which gives the result when a < 2(1 — fj.) (recall that w = \og{zg/z)). 

On the other hand, if a ^ 2(1 — ^), then the supremum in (2.70) is bounded for any value of 
w ^ 0, which shows the second part of the result. □ 



exponential convergence to equilibmum 25 

3. Spectral gap in weighted £^ spaces 
We now extend the above spectral gap result to the larger functional space X = (^{e^^Qi) given 

by 

X:={h^{h,), I :==^cxp(7?i)Q,|/i,| <+oo}, (3.1) 

■i 

for some < 77 < ^ log ^. The choice of rj ensures that this space is actually larger than the Hilbert 
space = £^{Q) used in the previous section, since 

/ 00 \ 1/2 / 00 \ 1/2 / 00 \ 1/2 

\\h\\x^ij2h^Q'] [J2Q,e^p{2r,i)j = V2 (j^Q^exp {2iji)j \\h\\^ , (3.2) 

and the latter parenthesis is finite for < ^ log ^. In order to see this notice that, since Qi = z^Qi, 

00 00 

Qi exp {2r]i) = Qj r' where r exp (2?; + log z) . 

i=l i=l 

Since the radius of convergence of the power series with coefficients i Qi is, by definition, Zg, the 
above sum is finite whenever r < Zs, this is, when ?7 < 5 log ^. 

Our approach to show that L can be extended to an unbounded operator on X that has a posi- 
tive spectral gap uses recent techniques [15] based upon a suitable decomposition of the linearized 
operator into a dissipative part and a "regularizing" part. We use the following result, which is a 
slight improvement over some of the consequences of [15]: 

Theorem 3.1. Let y, Z he two Banach spaces with norms \\ ■ \\y, \\ ■ \\z such that y ^ Z and 

\\h\\z^Cy\\h\\y, hey 

for some Cy > 0. Let Ly : 2>{Ly) -> y and Lz ■ ^(Lz) Z be unbounded operators in y and 
Z, respectively, and such that Lz is an extension of Ly (i.e., i^{Ly) C &{Lz) and Lz\si(Ly) ~ Ly). 
Assume that: 

(1) Ly generates a Co-semigroup {St)t^o on y. 

(2) It holds that 

WSthoWy ^Ci exp i-Xit) Wholly, hoey,t^O (3.3) 

for some Ci > 0, Ai e M. 

(3) Lz = A + B, where A, B are two linear operators on Z such that: 

(a) A: Z ^ y is bounded, i.e., 

IIA/ib ^ Ca||/i||2, h€Z, (3.4) 

for some Ca > 0. 

(b) B generates a Co-semigroup {S^)t^o on Z which satisfies 

llSfhollz^C2Cxp{-X2t)llhollz, hoeZ,t^O (3.5) 
for some C2 > 0, and some A2 > Ai. 
Then Lz generates a Co-semigroup (Vt)t^o on Z (an extension of {St)t^o) which satisfies 

llVthollz ^Cexp{-Xit)llhollz, hoeZ,t;^0 (3.6) 
for C = C2 + CyCiC2Cj^{\2 - 
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This theorem is directly inspired from the results in [15]. However, while the proofs in [15] are 
based on estimates of the resolvents of Ly and Lz, we present a completely different one based 
on the study of the dynamics generated by Ly and Lz- In addition to being remarkably simple, 
this proof allows us to state the result in general for any two Banach spaces (not necessarily one of 
them Hilbert) and to reach an exponential decay like Ai in (3.6). The resolvent methods in [15] 
are able to give more precise estimates on the behavior of the semigroup on Z, but here (3.6) will 
suffice. 

Proof. We notice first that A : Z Z is a bounded operator. Since B is the generator of a 
Co-semigroup (S'P)t^o in Z, by the bounded perturbation theorem, the operator Lz = A + B 
generates a Co-semigroup (Vt)tj>o in Z which, additionally, is given by the Duhamel's formula [4, 
Theorem 4.9] : 

VthQ = Sf ho+ f Vt-s (A Sfho) ds V/io e Z, t ^ 0. 



Thus, for fixed /iq e Z and t ^ 0, we have 

\\VM\z^\\Sfh4z+ I \\Vt-s{J^Sfh^)\\zds 
Jo 

^C2exp{~X2t)\\ho\\z+ [ \\Vt-s{ASfho)\\zds, (3.7) 

Jo 

where we used (3.5) in the last inequality. Now, for any s € (0, t), (A Sfho) G y since A maps Z 
to y. Then, since Lz is an extension of Ly, it is clear that (Vt)t^o is an extension of {St)t and 

Vt^s (A Sfho) = St-s (A Sfho) ey Vse (0, t). 
In particular, thanks to (3.3), 

/' WVt-s {ASfho) Wzds Cy f \\St-s{ASfh^)\\yds 

Jo Jo 

^CyCi f exp{-Xi{t-s))\\ASfho\\yds 
Jo 

^CyCiCA f exp{-Xi{t^s))\\Sfho\\zds. 
Jo 

Using again (3.5) we deduce from (3.7) 



. t 

-Alt I 



WVthoWz ^ C2Cxp{-X2t)\\ho\\z + CyC,C2C A e-^^' \\ho\\z / cxp(-(A2-Ai)s) ds 



^ C2Cxp(-A2t)||feo||z + ^y^'^f"^ cxp(-Ait) WhoWz 

M — Ai 



which, since A2 > Ai reads 



\yM\z ^ ( ^^'^'^f'^ + C2 I cxp(-Ai^) WhoWz, t^Q 



\ A2 — Ai 

and yields the desired result. □ 
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Our motivation for studying the spectral properties of the linearized operator in the larger space 
X is that in this larger space we have useful bounds of the nonlinear remainder term T{h, h): 

Proposition 3.2. Assume Hypothesis 2.1 and take < rj < ilog^. There is a constant C > 
depending only on {ai)i^i, and z such that 

\\r{h,h)\\x ^C\\h\\x\\h\\x, VheXi 

where Xi ~ i^iil + (Ji) exp{rii)Qi), i.e. 

Xi := {h^{h,), I \\h\\x, ■.^Y.^l + <J^)eMvi)Q^\h^\ < +c»} (3.8) 

i 

and we recall that {(Ji)i^i was defined in (2.16). 

Proof. Let us first treat the term Ti{h, h) (the first component oiV{h, h)). From (2.12) we have 

OO QO 

Qi\r,{h,h)\ = Qi\aiQihl + Y,a^Q^hlh^\ ^ aiQjhj + \hl\Y,a^QlQ^\h,\ 

i=l i=l 

CO 2 

.=1 ^1 

which bounds Ti{h, h). For the terms with i ^ 2, Ti{h, h) is given by (2.13): 

Ti{h,h) ^Tj{h,h) -Tl{h,h), with 

r.+ (/i, h) = fei = h, hi and TT{h,h) = a,Qih, (z > 2). 

We bound these terms separately. For r+(/i, h) = {T^{h, h))i^2 we have, using (1.5), 

OO CXD 

\\r+{h, h)\\x = J2 '^^P b,Q,\hi\ = \hi\Y, exp + 1)) h+iQ,+i\h,\ 
= exp(?7) \hi\ ^a,QiCxp(?7i)Q,|/i,| ^ cxp (77) \hi\ ^a^Q^cxp{?]i)\h,\ ^ ^^^^!2l\\h\\x \\h\\x,. 

The proof is even simpler for T~{h, h) = {T~ [h, h))i^2 since 

00 00 
\\r-{h,h)\\x = \hl\J2a^Qlexp{^Ji) Q,\h,\ < ^ a, exp (yyi) Q,\h,\ 

i=2 1=2 

which finishes the result. □ 



The following result shows that we may extend the linearized operator L from Definition 2.2 
to an operator on X with domain Xi. It can be proved by direct estimates on the expression of 
L given in (2.11) which are similar to those in the proof of the previous lemma, and we omit its 
proof: 
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Lemma 3.3. Assume Hypothesis 2.1 and (1.20c). There is a constant C > depending only on 
{ai)^^i, {bi)i^i and z such that 

\\L{h)\\x ^C\\h\\x, for all he £oo- 
Similarly, there is another constant C > depending only on {ai)i^i, {bi)i^i, z and r\ such that 

f\Q^^^^^^\Um<.C\\h\\x forallheioo- (3.9) 

The above lemma allows us to extend LtoXi: 

Definition 3.4. We denote by A the extension of the linearized operator L to the domain Xi. 

Theorem 3.5 (Extension of the spectral gap). Assume Hypothesis 2.1 and conditions (1.20a), 
(1.20c) and (1.20d), and assume also that z < Zg. Take < j] < ^log^. Then, the operator A 
generates a strongly continuous semigroup (exp(iA))t^o on X and there exists < A* ^ Aq such that, 
for some C > 0, 

oo 

|jcxp(tA)5||^ ^ C cxp{-X^,t)\\g\\x for allt^O and any g e X such that ^ iQigi = 0. 

i=l 

(We recall that \q is the size of the spectral gap of L in %, as given hy Definition 2.13.) In addition, if 
we assume 

lim Ui = +0O (3.10) 

i—^+oo 

then we may take A* = Aq. 

Proof We apply Theorem 3.1 to the Banach spaces 

oo 

■.= {hen\ ^^Q./l^ = 0}, 

1=1 

a subspace of "H, and 

oo 

x^ -.^{hex \J2iQ^h^ = 0}, 

a subspace of X. Note that under our conditions the map 

oo 

M : /ii — >M{h) ■.= ^iQihi 

i=l 

is continuous both in "H (as a consequence of (2.2)) and in X; hence, V.^ and X-^ are well-defined 
closed subspaces of "H and X, respectively. In analogy with (2.23) and (3.8) we define 

Hi ■.= H^nH2, Xi:=X^nXi. 

Let us briefly explain how one can project any sequence g £ X to X^. To do so, for any sequence 
g = (.g,),, let us introduce the sequence 

oo 

M^(g) = -M(g) Vi ^ 1 where V /cQ^. (3.11) 

^ fc=i 
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Since X]i>i cxp{r]i) Qi < oo, one checks that M : X ^ X isa bounded operator and for any g € X 
one easily checks that 

g-MgeX^. 

In the same way, the operator M : "H — > "H is bounded and h — Mh e "H^ for any h e %. 

To play the role of the linear operators Ly and Lz in Theorem 3.1 we take L-'- := L|^i with 
domain ^(L^) = and A-'- := K\x± with domain ^(A^) = X^. From (2.9) it is clear that the 
image of L is contained in Ti.^, and that of A is contained in X-^, so L-"- and A-"- are unbounded 
operators on "H^ and X^, respectively. 

Let us check the hypotheses of Theorem 3.1. As remarked in (3.2), %^ is contained in X^. 
Moreover, H-^ is dense in X^, since the set of compactly supported sequences with Y^^i iQih-i — 
is contained in both of them. Points (1) and (2) in Theorem 3.1 are established in Corollary 2.10 
where it is shown that one may take Ci = 1 and Ai = Ao in Theorem 3.1, where Aq is the spectral 
gap estimated in Section 2.3. 

It remains to show point (3) of Theorem 3.1. 

Step 1: the splitting A^ = A + B. Take an integer iV ^ 1, to be fixed later. We first split A as 
A + B'ma. similar way to Proposition 2.9, and we then write A-^ = + by projecting to X^. 
We define L'^ and L'^ similarly to (2.31)-(2.32), but this time in X: we choose iV > 2 and set, for 
all i ^ 1, 

OO 

(/i) := '^X{uAn{i,]}^N} ^ijhj - cJihiX{i<^N} for heX, (3.12) 

OO 

^^(/i) -o^hiX{t>N} + X] X{min{i.j]>N} £.i,]hj for h e Xi. (3.13) 

We will see below that the sums converge absolutely for h in the corresponding domain, so the 
definition is correct. We also define the truncated identity operator = as 

and then set, for i? > to be fixed later, 

A^h!^ + RI^, hex 
B = A^'^-RI^, heXi, 

so that A = A + B. We finally project this splitting to X^ thanks to the operator M introduced in 
(3.11). Namely let 

A^ h = Ah- M{Ah), h e X^ 

B^h = Bh-M{Bh), heX^, 

Step 2: A is "regularizing". We show now that Ah eH for all h e X (which, since M : H ^ H 

is bounded will imply then that A^{h) e "H^ for all h e X^). Notice that one may write, more 
explicitly, 

^iih)=j:jL2^i,jhj-aih, A^{h)^^2Ah, + a2Qih3~a2h2 
^?{h) = £.i,ihi + hi-i X{iiiN+i} + a.iQi /i^+i X{i^w} - cfihiX{iiC,N] for i > 2. 
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The only infinite sum in the definition of A (the one in Af ) converges since, from (2.17)-(2.18), 

oo ^ oo oo 2 

In order to see that A{h) e H it is enough to show that J2i1n+2 Qi-^iW^ =^ for some 

C > (we omit the first + 1 terms, for which a similar bound is obviously true) . But this is true 
since for i ;^ + 2 the only term in Ai{h) is the first one in the expression of Af{h) in (3.14) and 

oo oo oo . 

i=N+2 i=N+2 i=N+2 1 

where A is defined in (2.2). Thus, A : X ^ H is bounded, and hence A^ : also is. 

Step 3: S is strictly dissipative. We prove that 

(sign(/i),6/i)^, ^ < -X3\\h\\x for all h e Xi, (3.15) 

where (•, ■)x',x denotes the duality pairing between X and its dual X', and A3 will be any number 
larger than Ao when (3.10) holds, and some positive number otherwise. One sees that (3.13) may 
be rewritten as 

Af^{h) = ~(JihiX{i>N} + bthi^iX{i>N+i} + atQih,+iX{i>N} (« ^ !)• (3.16) 
For any h & Xi one has 

00 00 AT 

{sign{h),Bh)-^,^^ = ^ cxp (771) Qi sign{hi)Bi{h) = - ^ cxp criQi\h,\-R'^exp Qi\h 

i=l i=N+l i=l 

00 00 

+ E "sxp (771) Qi sign{hi) fej/ii-i + E "^^P (^*) QiQi sign(/ij) a.hi+i. 

i=N+l i=N 

The last two sums can be bounded by 

00 oo 

E exp (jji) Qibi\hi-i\ + E exp irji) QiQiai\hi+i\ 

i=N+l i=N 
7i Qi 



^ ( cxp(?7) —^h+i +cxp(-?7) ^ ' ^ Q»-i ) cxp (771)0-^2^1/1^1 = ^iiCx:p{rii)aiQi\K\ 

i=N ^ ^ z— iV 

where 

/ij := cxp (7/) — !— ^ + cxp (-7/) — . (3.17) 

Now, due to (1.20c) and (1.20d), has a limit as i +cx): 

Qia, Qia^ Oi 1 



where I := Zg/z ;^ 1, while 



EXPONENTIAL CONVERGENCE TO EQUILIBMUM 



31 



Hence, 



exp (rj) + cxp {—1]) £ 



as z +00. 



One easily checks that this hmit is strictly less than 1 if and only if exp (77) < £. Since we are 
assuming exp (77) < \/I ^ i, for N large enough we can find /i < 1 such that 



N 



(sign(/i),i5/i)^,_^ ^ -(1 - /^) ^ exp(?7i)CTjQi|/i,| - i?^exp(?7i) Qi|/ii 

i=N+l i=l 

Choosing A3 = (1 — ^,)a and i? = A3 we obtain 



{sign{h),Bh)x, x < -A3 ^ cxp(?7i)Qi|/ii 



This is nothing but (3.15). 

If we assume additionally that (3.10) holds (i.e., linij^oo clj ~ +00) then take any A3 > Aq. 
One can choose iV > 1 large enough so that ct^ > for any i ^ N and then, by picking i? ^ A3 
we get 

00 



1=1 



which proves (3.15) with any A3 > Aq. 



Step 4: is strictly dissipative. We now show that if we take any A2 < A3 (with the A3 from 
Step 3) then one can additionally choose N in (3.12)-(3.13) such that 



{sign{h),B^h) ,.,_^ < -X2\\h\\x for all h e X^. 



(3.18) 



This is not obvious, since the property (3.15) of B is not necessarily shared by its projection B^. 
However, in this case and for N large enough, B-^ happens to be a small perturbation of B on X^. 
We need to estimate, for h £ X^, the quantity 



N 



m{Bh)^ J2 ^Q^^^(h)-RY,lQ,K 

i=N+l 



i=l 



We have, with very similar calculations as those needed to obtain (3.9), that 

00 OC / .\ 

E ^Q^\^^f{h)\^e,{N) Q,!^£^|Af(/^)KCel(7V)||/^|U 



(3.19) 



with 



i=N+l 



ei{N) := sup — -, 

i^N+i exp(77i) 



which tends to as — > +00 (as implied by (1.20c) and (1.20d) by considering the ratio of two 
consecutive terms in the sequence). On the other hand, since J2'^i iQi^i = 0, 



R 



N 



R 



i=N+l 



^ -ei{N)\\h\\x 
a 



(3.20) 
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where we recall that = inf j^i o-^. From (3.19), (3.20), (3.11) and (3.15) we have 

{sign{h),B^h)^,^_^ < _A3||/.||x + Ei2i£^pM +c) e,{N)\\h\U 

for all h e X^. Choosing N large enough proves (3.18), since ei(A^) -> as -> +oo (notice that 
the choice of R from Step 3 is not affected.) 

Step 5: B generates a semigroup on X. To prove that B is the generator of a Co-semigroup on 
X we shall invoke Miyadera perturbation Theorem [4] . Notice that B may be written as 

B^T + C 

where T is the multiplication operator (with domain S'{T) = S>{B) = Xi) given by 

%{h) = - (a, + R) X{^^N}h^ Wi^l, he S!{B) 

while C is defined by 

Ct{h) = hhi^iX{i>N+i} + at Qi ht+iX{i>N}, « > 1- 

It is clear that T is the generator of a Co-semigroup (U*)f^o on X given by 

Vl{h) =cxp{-t{a^ + R)x{t;,N}) i ^ 1, t ^ 0, heX. 

Moreover, C : Xi ^ X is bounded. Indeed, it is not difficult to check, as we did in the proof of 
(3.15) that 

OO 

\\Cih)\\x ^ fJ'iCTiexpiiji) Qi\hi\ ^ \\h\\xi sup /i^ \fh ^ {hi)i^i e Xi (3.21) 

where is defined by (3.17). Since the sequence is bounded, this shows in particular, that C 
is T-bounded. One deduces then from (3.21) that 

oo 

||CU*/ii|A' s$ X! /^^c^^exp(-^(a, + i?))exp(77i)Q,|/i,| h e X^ , t ^ 

i=N+l 

which readily yields 

/ \\CV'h\\xdt^ V L^^expirji) Q,\h,\. 

In particular, we see that 

/ \\C\5*h\\xdti^^i\\h\\x yheX^ 
Jo 

where yu := sup^^x+i l-^i- already saw in Step 3 that N can be chosen large enough so that 
^ < 1. Then, the above inequality exactly means that C is a Miyadera perturbation of T (see [4, 
Section 4.4, p. 127-128]) and that S = T + C is the generator of a Co-semigroup (S'f )t^o on X. 
Notice then that (3.15) exactly means that (A3 + B) is the generator of a contraction semigroups 
inX. 

Step 6: (A2 + B-^) generates a Co-semigroup of contractions in X-^ Notice that the approach 
used in the previous step seems difficult to apply directly to B^, the reason being essentially that, 
in the above splitting B = T + C, the operators T and C does not map their respective domains to 
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X^. To prove that (A2 + B-^) generates a Co-semigroup of contractions in X^, we adopt another 
more indirect v^fay. Since A : X ^ Ti is bounded, it is clear that A : X X is a bounded 
operator. Therefore, by the bounded perturbation theorem, A = A + B is the generator of a Co 
semigroup {Tt)t^o in X. Moreover, since A{Xi) c X-^, the closed subspace X^ c X is invariant 
under {Tt)t^o- There, as well-known [10, Chapter II. 2. 3, p. 60-61], the restriction A-'- = A\x± 
with domain &{A-^) ~ Si(A) n X^ = X^ is the generator of a Co-semigroup in X^. We already 
saw that A-^ = A^ + B-^ where B^ is a bounded operator with domain X^ and A-^ : X^ 
is bounded. In particular, A-^ : X ^ X is also a bounded operator and then, 

B^ = -A^ + A-L 

is the bounded perturbation of the generator A^. Thus, B-^ is the generator of a Co -semigroup 
[Sf )t in X^. Moreover, according to Step 4, inequality (3.18) implies that A2 + is a dissi- 
pative operator in X^ (in the sense of [10, Chapter II.3.b, p. 82]). Since A2 + B^ generates a 
Co-semigroup in X-^, according to the Lumer-Phillips Theorem, this semigroup is a contraction 
semigroup, or equivalently 

WsfhaWx i^e^v{-^2t)\\ho\\x ^hoeX^, t^O 
i.e. (3.5) holds with C2 = 1 and A2 provided by (3.18). 

Step 7: Conclusion. To obtain the desired conclusion, we only have to apply Theorem 3.1 with 
Ai = Ao being the spectral gap of L while < A* = A2 is any positive number strictly smaller than 
min(Ao, A3) where A3 is constructed in Step 3. We just recall here that, if limi^oo ai = +00, the 
number A3 is any arbitrary real number strictly larger than Aq, and thus in this case any A* S (0, Aq) 
would yield the conclusion. □ 

4. Exponential convergence for the Becker-Doring equations 

4.1. Local exponential convergence. In this section we prove a local version of Theorem 1.1. 

Theorem 4.1 (Local exponential convergence for the Becker-Doring equations). Assume (1.19) 
and (1.20). Let c = {ci)i be a nonzero subcritical solution to equation (1.1) (i.e., with density g < gj 
with initial condition c(0) such that (1.22) holds for some ly > 0. Take z > satisjying (1.7). 
Then there exist some < 7] < i^, some C, e > and some A* > (all explicit) such that if 

00 

^exp(?7i) \c^{0) ~Qi\<e 

i=l 

then 

oo 

^cxp(ryi)|Q(t) - Q.,| s; Ccxp(-A,i) forallf^O. (4.1) 

1=1 

In addition, A* can be taken equal to Ao (the size of the spectral gap of the linear operator L) if 
a-t = +00. 

We remark that all the constants in the above result can be explicitly estimated. The proof of 
the above local convergence results relies on two crucial estimates. The first one was proved in 
Prop. 3.2, a bound of the nonlinear term T{h, h); the second one is a uniform-in-time bound of 
exponential moments for the Becker-Doring equations, which we take from [16] : 
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Proposition 4.2. Assume (1.19) and the conditions (1.20). Let c = (c,),; be a nonzero subcritical 
solution to equation (1.1) with initial condition c(0) such that (1.22) holds for some v > 0. Then 
there exists < F < and Ki > such that 

oo 

^exp(I7j)cj (t) A'l, (4.2) 

where c = {ci{t))i is the unique solution to (1.1) with initial datum c(0) = (ci(0))i^i. 

We notice that (1.19) and (1.20), together with the definition of Qi through (1.5), imply condi- 
tions (H1)-(H4) of [16]. In particular, since Zg is defined as the radius of convergence of the series 
z^iQi, it is clear that 

lim Q]'' = 

which is (H3) of [16]. 

We are ready then to prove Theorem 4.1: 

Proof of Theorem 4.1. Let c be a solution of (1.1) and let h = {hj) be the fluctuation around the 
equilibrium, defined as in (2.4). Since Proposition 4.2 holds under our conditions, we can find 
Q < V < v such that (4.2) holds. We take any 7] < min{77, i log ^}, and consider X the vector 
space defined in (3.1). Note that Proposition 3.2 and Theorem 3.5 are applicable under these 
conditions, and consider the quantity appearing in Theorem 3.5. 
The fluctuation h satisfies equation (2.8) in X: 

-^h^k[h]+T{h,h), 

so, if {St)t^o denotes the evolution semigroup generated by A, then 

h{t)^Sth{Q)+ f St-sT{h{s),h{s))ds. 
Jo 

Recall that X]i>i ihi{t)Qi = = X^i^i i^i{h{t), h{t))Qi for any t ^ so that, according to Proposi- 
tion 3.2 and Theorem 3.5, we have (for some constant C > 0), 

\\h{t)\\x ^\\Sth{0)\\x+ f \\St-sr(h(s),h{s))\\xds 
Jo 

=^ C\\h{0)\\x cxp(-A*i) + Ccxp(-A^t) / cxp(A*s)||r(/i(s), h{s))\\xds 



^C\\h{0)\\xexp{-\,t) + Ccxp{-\J) / expiX,s)\\his)\\x\\his)\\x^ds (4.3) 

Jo 

where we recall that Xa is defined by (3.8). For any 6 e (0, 77), we have from Cauchy-Schwarz's 
inequality that 

IIMi)llA-„ =^|/j,(i)N"Q.exp(7;*) 

i 

/ 00 \ 1/2 / 00 \ 1/2 

^ K]|/^,(^)|^"Q^exp(^?^-<5^) ^ |/i,(ON"Q. exp (771 -f <5z) . 
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Now, using that 

i" exp (—Si) ^ C, i" exp {r]i + Si) ^ C cxp (Vi) for all i ^ 1 

we deduce that 

/ oo s 1/2 

\\h{t)\\x^^C\\h{t)\\][^(J2Mt)\Q,cxp{m)\ ^CKl^^\\h{t)\\]l^ foralH^O. 
^ 1=1 ^ 

Using this in eq. (4.3) we have 

pt 

\\h{t)\\x ^C\\h{0)\\xexp{-X,t) + CeM-^*t) / exp(A^s)||/i(s)||^/' ds. 

Jo 

Then, the quantity := exp{X^,t)\\h{t)\\x satisfies 

u{t)i:C\\h{0)\\x + C f ii(s)3/2ds. 

Jq 

We deduce then the conclusion from Gronwall's inequality. □ 

4.2. Global exponential convergence. In order to prove Theorem 1.1 we complement the local 
result given in Theorem 4.1 with the following global result from [16]: 

Proposition 4.3. Assume (1.19) and the conditions (1.20). Let c = {ci)i be a nonzero subcritical 
solution to equation (1.1) with initial condition c(0) such that (1.22) holds for some v > 0, and take 
z > satisfying (1.7) as usual. 

Then, there exist two explicit constants Ci > and ki > such that 

|cj W - Qj K Ci exp(-Kii^) Vi ^ 0. (4.4) 

Proof of Theorem 1.1. Let h = {hi) be defined as in (2.4). Let us fix the same r/ e (0,17) as in 
Theorem 4.1 and denote by X the £^ space with weight e^'* that was defined in (3.1). The idea is 
to use Proposition 4.3 to show that at a certain time Iq we are close enough to the equilibrium to 
use Theorem 4.1. To do this we use a simple interpolation argument to estimate the norm 

Y,\c^{t)-Q,\exp{T^i)^\\h{t)\\x 

in terms of X]j>i * {t) — Qi I and a slightly stronger norm. Namely, for any S € (0, 77), a simple 
interpolation argument yields 



i^i 



< ^i|/i.(t)|Q, K]|/i,(t)|Q,cxpC 
Therefore, according to (4.2) and (4.4), one obtains 

S + T] 



hit)\\x ^ C/+" cxp ( -^^ti ] yt^O. (4.5) 
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Consider the number e > given by Theorem 4.1. Then, there exists an exphcit to > such that 

II'^(^o)ILy s$ e 

and one deduces from (4.1) starting from t = to that 

\\h{t)\\x ^ CJh{to)\\x exp(-A,t) V< ^ to- 
Together with (4.5) for t < to, this concludes the proof. □ 

Appendix A. Discrete Hardy's inequalities 

The following discrete version of Hardy's inequality is used in the proof of Theorem 2.14: 

Theorem A. 1. Consider two sequences of positive numbers {pLi)i and (i^i);. Then, the following are 
equivalent: 

(1 ) There exists a finite constant such that 

oo / i \ oo 

i=l / i=l 

for any sequence f = (fi),. 

(2) The following holds 



< oo. 



If these equivalent propositions hold true, then B A ^ AB. 

Remark A. 2. The above discrete version of Hardy's inequality is not new [19] and we give a proof 
for the sake of completeness. It is a simple adaptation of the proof of the continuous version of Hardy's 
inequality, originally due to Muckenhoupt [21], that can be found in [1, Chapitre 6]. 

Proof Assume that there exists some finite A > such that (A.l) hold true for any / = {fi)i and 
let us prove that A^ B. For any fixed k ^ 1, set 

{— for i ^ fc 
for i > k. 
Then, one recognizes easily that X^i^i '^iff = J2i=i T while 




i=l 

According to (A.1) one gets then 

k / oo \ / k 



i=l ' \i=k / \j=l ^ 
or equivalently A > (X^i^fc V^i) 17") f^'^ fc ^ 1. This proves that B ^ A. 
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Conversely, let us assume that B < oo and let us prove that A ^ 45. Set 

Ti = ^ — , T, = ^nj, 7j = and /3j = \/T~. 

Let / = {fi)i be a given sequence. Thanks to Cauchy-Schv^arz inequality 

2 



Now, since ^ = — Tj_i, one gets 

1 1 _ T, - T,_i 



Using the general inequality 

X -Y 



we get that 



s; 2 f - V7 ) yx >o 



Therefore, the left-hand-side of (A.1) can be estimated as follows: 

oo/iX oo / \ ^ 

^2^/i,A E/'^^/^O =2E/>^/3, |Em./3. 

i=l = l / 1=1 \3 = 1 / j=l 

Then, one sees that to prove our claim, it suffices to prove that 



One notices that, by definition of _B, /3i^i < vi3 for any i ^ 1. Therefore, 

oo oo oo 

fi ■ - V -L i 

i=j j=j i=j 

Now, since — F^+i — fii one has, from (A.2): 



and the above inequality reads 



^ 2 Vrl- v/rw 



E M« ft < 2/bE ( v/tI - VrvT) = 2^/5^17 = Vj ^ 1. 
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(B.l) 



Finally, by definition of B, we get that 

which gives (A.3) and achieves the proof. □ 

Appendix B. Some estimates of special functions 

We collect here some technical estimates for special functions which turn useful for estimating 
the spectral gap size near the critical density (see Section 2.4). 

Lemma B. 1. Take < ^ ^ 1. For some quantity depending only on ^, 

/>oo 

/ expi-yf") dy ^ exp{-x'')x^^f' for all x ^ 1 

J X 

Proof. Through the change of variables u = y^^ — x'^ we obtain 

f°° 1 f°° / u \ ^ 

/ exp(-y^) dy = — cxp(— a:;^)a;^^'^ / cxp(— u) I hi) ^ du 

Jx Jo ^a;/^ ^ 

so the result holds with = ^ /g°° cxp(— m) (u + du. □ 

Lemma B. 2. Take < /i ^ 1 and < a < 1. For some quantity C^,a depending only on fi and a, 

[ cxp{y^')y-°' dy C^,^ cxp{x^')x^~^'-" for all x ^ 1. (B.2) 
Jo 

Proof We first break the integral in the intervals (0, x/2] and {x/2, x): for the first part, 

exp(y'^)y-"dy cxp((x-/2)^) ^ ' y-" dy = -1- exp((a;/2)^) (|)'"" ^ cxp(x^)a; 



(B.3) 



with 

2a — 1 

Af,,a z max{exp((a;/2)^ - . 

1 — a x^i 

For the second part, through the change of variables u = — we obtain 



exp(y^)y-" = ^ cxp{-x^')x^-^'-°' ' exp(-u) - ^) 



du 



^ exp(-a;^)a;i-^-" / ' cxp{~u) du B^,,aexp{-x^')x^-^'-°' (B.4) 

Jo 



with 

:= -max|l,2i-^^-° | 

since |<l-^^lon the region of integration. From (B.3) and (B.4) we conclude with 
We end this Appendix with an asymptotic result for special expansion series: 
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Lemma B. 3. Take fi > and a sequence {gk)k such that, for some Cq > 0, 

gk+i - gk ^ Co/s''"^ ask^ +00. (B.5) 

Then 

k 

y^eMgjW ~ TT exp{gk)k''+'-^ as k ^ +oo (B.6) 

and 

oo ^ 

V exp(-5,)j'' ^ — exp(-5fc)fc'+'"'' as k ^ +cx3. (B.7) 

j = k 

Remark B. 4. One easily deduces that, in the conditions of the above lemma, there is a constant C 
which depends only on the sequence {gk)k arid v such that 

k oo 

^exp(g,)/ ^Ccxp(gfc)fc''+i-^ and ^ exp(-g,)r Cexp(-gfc)fc'^+'"'' 

j=l j=k 

for all k^l. 

Proof We prove (B.6) by using the Stolz-Cesaro theorem (which is a discrete analogue of THopital's 
rule). Namely, setting 

k 

iJfe = exp(gfc)r+i-'' and Vk = J2exp{gj)f (fc ^ 1), 
in order to show that 

hm — = Go, 

k-foo Vk 

lim = Co. (B.8) 

fc^oo Vk+i — Vk 



we instead prove that 



One has 

Uk+i - Uk ^ (k + - cxp(gfc - fffc+i)fc-+i-^ 

Vk+i - 14 (fc + 1)^ 

_ ((fc + 1)"+!-^' - + fc-^+i-^' (1 ~ exp(,gfc - fffc+i)) 

(fc + 1)'' 

Since 

|(fc + l)'-+i-A--fc-+i-A'| |t. + l-^|min{(fc+l)'-^fc'-^} 

^ 7-. ^ + 



{k + iy {k + iY 

as /c ^ +00 while 

fc-+i-^(l-cxp(gfc-fffc+i)) fc^+i-A-Cgfe+i - fffe) 1 - cxp(gfc - gk+i) „ 
7- — = lim = Go, 

k^oo + fc^oo (fc + 1)'' gk+i-gk 

using (1 — exp(gfe — gk+i)) /(.9fe+i — gk) -> 1 together with (B.5), we obtain (B.8). We apply the 
Stolz-Cesaro theorem and obtain (B.6) (notice that, since (B.5) implies that gk ^ (G'o//i)fc^, the 
sum J2'j=i cxp(gj)j'' diverges and hence the hypotheses of the theorem are satisfied). Equation 
(B.7) can be proved by a very similar calculation. □ 
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